Non-Fermi-Liquid Scaling in Ce(Ru .5Rho.5)2SÌ2 
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We study the temperature and field dependence of the magnetic and transport properties of the 
non-Fermi-liquid compound Ce(Ruo.5Rho.5)2SÌ2- For fields H < 0.1 T the results suggest that the 
observed NFL behavior is disorder-driven. For higher fields, however, magnetic and transport prop- 
erties are dominated by the coupling of the conduction electrons to criticai spin-fiuctuations. The 
temperature dependence of the susceptibility as well as the scaling properties of the magnetoresis- 
tance are in very good agreement with the predictions of recent dynamical mean-field theories of 
Kondo alloys close to a spin-glass quantum criticai point. 
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The properties of a large class of /-electrons materials 
show striking departures from the predictions of stan- 
dard Fermi- liquid theory at low temperature Q] . Several 
mechanisms leading to non-Fermi-liquid (NFL) behavior 
have been proposed. In systems close to a quantum phase 
transition such as CeCu^ x A\i x M or Celn 3 ||, NFL be- 
havior is due to the coupling of the conduction electrons 
to criticai spin fluctuations Q . Anomalous properties are 
observed when the system is driven through the quantum 
criticai point (QCP) by alloying or by applying pressure. 
In other systems, such as UCus-^Pd.,; ||, NFL proper- 
ties are thought to be a consequence of the interplay of 
strong structural disorder and many body effects 
In this paper we report results of a study of the temper- 
ature and field dependence of the magnetic and trans- 
port properties of Ce(Ruo.5Rho.s)2SÌ2. We found that 
this system exhibits different types of anomalies depend- 
ing on the value of the applied field. At weak fields, 
we found signatures of disorder effects such as a diverg- 
ing low-temperature susceptibility and an anomaly in the 
low-field magnetoresistance. Above lkG, the T and H- 
dependence of the susceptibility and the resistivity agrees 
with the predictions of recent mean-field theories of the 
spin-glass (SG) QCP The magnetoresistance is 

found to exhibit universal scaling properties as predicted 
by the theory. 

In the Ce(Rui_ a; Rh a ;)2SÌ2 alloy series the Ce-sublattice 
is preserved and the hybridization between 4/ and con- 
duction electrons varies with the concentration of the 
ligand Ad atoms Ru and Rh. Pure CeRu2SÌ2 is a 
heavy fermion compound with a 7 value of about 385 
mJ/mol/K 2 |Tl"| and no long range magnetic order down 
to 20 mK. With substitution of Rh for Ru, a spin den- 
sity wave (SDW) region appears between x — 0.03 and 
x = 0.4 @[n§. While Ce(Rui_ 2; Rh K )2SÌ2 at x = 0.03 
is a normal Fermi liquid [jl4 jl5 1 , a NFL regime exists for 
x = 0.4 and 0.5 lHQ. In pure CeRh 2 Si 2 , the 4/ elec- 
trons are localized and the material is antiferromagnetic 
(AF) below T N = 35 K |l7j]. With increasing Ru substi- 
tution, TV decreases and eventually vanishes at a criticai 



concentration x c ~ 0.5 |18| . There is no direct evidence of 
long range AF order below x = 0.7. The admixture of Ru 
and Rh introduces magnetic frustration effects as it leads 
to a compctition between widely different types of mag- 
netic short-range order ]j~6| ]. A frustrated ground state 
of the SG type can not be excluded slightly above x c . 
Recent ^-SR studies [|r| showed that the T-dependence 
of the muon relaxation rate in the x = 0.5 alloy is similar 
to that observed in spin glasses. The muon depolariza- 
tion rate decays exponentially, however, showing that the 
spin correlations are dynamic rather than static. 

Samples of Ce(Ruo.5Rho.s)2SÌ2 were prepared by arc- 
melting of the constituents in an argon atmosphere. The 
ingots were remelted several times to insure homogene- 
ity. Single crystals oriented along a and c-axis were 
grown by the Czochralsky method in a tri-are furnace 
in argon atmosphere and parallelepiped-shaped samples 
of size «0. 5x0. 5x4 mm 3 were obtained. The resistiv- 
ity was measured by a standard ac method at 17 Hz in 
the range 16 mK - 4 K and in magnetic fields up to 5 T. 
The low-field susceptibility was measured in the dilution 
range ss 50 mK - 4 K using a standard method at 130 
Hz. The static magnetization was measured in the same 
range of temperature with a SQUID magnetometer in a 
dilution setup and, above 3 K, in a commercial SQUID 
magnetometer . 

Fig. [l] displays the resistivity p measured in a field 
applied along the c-direction up to 5 T. The data corre- 
spond to a current flow in the a-direction. The behavior 
along the c-axis is similar but the resistivity is about 
four times smaller. The magnetoresistance is positive at 
low temperatures and changes sign at about 2.5 K in 
the a-direction and at about 1.5 K in the c-direction. In 
both cases the data at 5 T follow a T 2 -law up to T rj 
1 K. The range of temperatures in which FL behavior 
is observed decreases with field and vanishes at H = 0. 
These results are qualitatively similar to those obtained 
in CeCu6~xAu x J2J and in the stoichiometric compound 
CeNÌ2Ge2 J2(|. The resistivity of our samples is much 
higher, however, as a consequence of the high degree of 



1 



substitutional disorder present in the x — 0.5 alloy. 
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FIG. 1. Resistivity of Ce(Ruo.5Rho. 5)2812 in a magnetic 
field as a function of T 2 . The magnetic field is along the 
c-axis and the current flows along the a-axis. The arrows 
indicate the range of temperature where the resistivity varies 
as T 2 . Inset: field-dependence of A = dp/dT 2 determined 
from the data below 100 mK. 

The temperature dependence of the zero-field resistiv- 
ity is 5p oc T 16 as shown in Fig. ||(a). The resistivity 
exponent is close but not identical (see below) to that ex- 
pected for metallic antiferromagnets Q and spin-glasses 
[§,§ at the QCP in the high resistivity limit [gì}. The 
low-temperature magnetoresistance (p(H) — p(0))/p(Q) 
along the a-direction is plotted in Fig. ||(b) for H < 1 T. 
At high fields (not shown in the figure) it shows the clas- 
sical ff 2 -dependence due to the bending of the electron 
orbits. Below a few kG, however, the low-temperature 
resistance increases linearly with field as H — > 0. 

The susceptibility x = M/H (M is the magnetization) 
is represented in Fig. |3| for H ac = 1 G and H^ c = 0.01 
T, 0.1 T, and 1 T for T < 10 K. At low temperatures \ 
decreases strongly with increasing field between 1 G and 
1 kG but weakly above 1 kG. Above 3 K the field has no 
effect on x up to 1 T as seen by comparison of the curves 
at 1 kG and 1 T in Fig. (|3|). At 1 G, \ increases sharply 
with decreasing T below 2 K. 

The divergence of \ a t low-fields as well as that of 
C/T [H| suggests that, in this regime, NFL behavior 
may be driven by disorder j6[0] . The linear (rather than 
quadratic) rise of the low-field magnetoresistance at low 
temperatures reported here may be explained by Kondo- 
disorder effects ^] . Recent /zSR experiments |23| showed 
a sharp increase in the muon relaxation rate below 2 



K that saturates at a T-independent value below 0.7 
K. This behavior has been interpreted in terms of the 
Griffiths-phase scenario in which finite clusters carry- 
ing a magnetic moment fluctuate very slowly at low tem- 
perature due to quantum tunneling. Magnetic, NMR 
and specific heat experiments performed in a temper- 
autre range much higher than ours have also been inter- 
preted in terms of this mechanism H.p5| . 
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FIG. 2. (a) Temperature dependence of the zero-field resis- 
tivity of Ce(Ruo.5Rho. 5)2812 along the a-direction. (b) Magne- 
toresistance p(H) — p(0))/p(0) measured along the a-direction 
at several temperatures and plotted as a function of the field 
applied in the c-direction. 

For H > 1 kG, \ remains finite as T — > and depends 
weakly on H, suggesting that application of a moder- 
ate magnetic field quenches the mechanism that leads 
to the divergence observed at lower fields. Although a 
full understanding of this fact is stili lacking, it should 
be noticed that in the Griffiths-phase model M quan- 
tum fluctuations of the largest clusters are expected to 
be suppressed by a small magnetic field. Indeed, while 
the Zeeman energy of a cluster of size N grows as vJV, 
its tunneling energy vanishes cxponentially with N. In 
the following we concentrate on the physics above 1 kG. 

The T-dependence of the susceptibility at 1 kG is stili 
anomalous and Sx(T) is approximately linear in T 3 / 4 (cf. 
Fig. ||) except at the lowest temperatures, a point that 
will be discussed further below ^6|. A T 3 / 4 dependence 
of x as wei l as the value 3/2 for the resistivity exponent 
where predicted by recent dynamical mean field-theories 
(DMFT) of the spin-glass QCP |-[l0|- In vicw of the im- 
portami role that frustration is expected to play in this 
compound it is tempting to try to interpret our results 
in terms of the fully frustrated SG model. The latter 
describes conduction band electrons coupled to localized 
/-electron spins via a locai Kondo coupling, Jk- Thcrc 
is, in addition, a residuai Ising-like interaction between 
the localized spins. The effects of the magnetic frustra- 
tion introduced by disorder are incorporated by taking 
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random spin couplings Jìj chosen from a symmetric dis- 
tribution of width (Jfj) = J 2 ■ Disorder in the Kondo 
temperature is not included in the model. From the 
Kondo-temperature distribution determined in Refs. |24| ] 
and [^|| one can conclude that this effect should play a 
lesser role than frustration in the low-temperature range 
that interests us. The SG model |)|,[l0) was investigated in 
the framework of dynamical mean-field theory (27| . It ex- 
hibits a zero-temperature SG transition when the typical 
exchangc encrgy J = J c ~ Tk, the Kondo temperature of 
the underlying Kondo lattice. Monte Carlo simulations 
pp| | of this model showed that its criticai properties are 
described by an effective strong-coupling theory closely 
related to other mean-field models f§|J|. The physical 
properties of the system depend on the effective distance 
to the QCP, A(T, H). It can be shown that this is 
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where A = (1 - J/ J c ) + {H/H f and H = J c /(gm) 
(g is the gyromagnetic ratio of the Ce ion). The scale 
To is proportional to Tk- Numerical simulations yield 
J c ~ 1.15 Tk- The spin susceptibility is pi] 



J C X = VT+ A - VA. 



(2) 



The spectrum of magnetic excitation has a scaling form 
JcX"(u) = %/A$(w/J c A), where the universal 
scaling function $(x) = x/V2 [(1 + x 2 ) 1 ' 2 + l]" 1 / 2 . The 
temperature-dependent contribution to the resistivity is 
computed from Sp oc 1/r with the inverse scattering time 
r -1 ex J Q °° <Ajx"(w)/ s'mh(f3uj), an expression valid in the 
dirty limit plj . The resistivity has the scaling form 



p(T,H)- p(0,H) cxT 3 / 2 \f 
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with = x x / 2 du<& (ux) / sinhu. It follows that 

Sx = X(0)-X(T) oc T 3 / 4 and Sp oc T 3 / 2 at the QCP. The 
resistivity exponent of the mean-field SG model coincides 
with that of the d = 3, z = 2 antiferromagnet ||||]. The 
susceptibility exponent is specific to the SG model. Away 
from the QCP, normal Fermi-liquid behavior as T — > is 
recovered with both Sx and Sp oc T 2 /\/Aq at sufficiently 
low T. The crossover between these limiting forms will 
be discussed below. The parameters of the theory can be 
determined from an analysis of the experimental data. At 
the criticai concentration, r = 1 — J j ' J c vanishes. How- 
ever, x c is not known accurately and a small but finite r 
can not be excluded a priori. The characteristic field Hq 
may be estimated from an extrapolation to T = of the 
susceptibility per Ce atom. From the definitions above 
and Eq. (|^) we estimate H = /ìb/x(0) ~ 11 T (we have 
assumed g = 2). Since Hq is very large, the measuring 
field can be neglected in the analysis of the magnetiza- 
tion at 1 kG. A fit of x using Eqs. (||) and (|l|) with 



Aq = gives To « 20K, which is slightly smaller than 
the Kondo temperature of the system estimated from the 
T-dependence of the resistivity |15[ . 
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FIG. 3. Susceptibility M/H measured in several magnetic 
fields along the c direction, plotted as a function of T 3/4 . The 
solid curve is a fit of the data at H=l kG to the expressions 
in Eqs. (^) and (^). Inset: The ac susceptibility at H=l G. 
The solid line is the fit mentioned in the text. 

We analyzed the T- and if-dependence of the resistiv- 
ity using Eqs. (||) and ([j]). The condition that ali the data 
in Fig. ^] collapse into a single scaling curve leaves little 
freedom in the choice of the parameters. In particular, 
we found it impossible to scale ali the curves with r = 0. 
A scaling plot of the resistivity along the a-axis is shown 
in Fig. [i[ The data points are the values of the scaled 
resistance (p(T, H) - p(0, H)) x T~ 3 / 2 plotted vs the re- 
duced variable t/A(t = T/T ) for T < 0.9 K and H < 5 
T. The values of the parameters are r = 7 x 10~ 3 , T = 20 
K and H = 13 T. The value of r measures the distance 
to the true QCP, r = 5J/J C giving 5J « 0.2 K, a very 
small energy compared to the other energy scales present 
in the problem. The characteristic field determined from 
this analysis is close to the theoretical estimate given 
above. The solid line in Fig. ^ is the theoretical scal- 
ing function ^>(x). There are no adjustable parameters 
other than an amplitude that fixes the vertical scale. The 
agreement between theory and experiment is very good 
except for the data for H=0 (the empty squares in Fig. 
||) which lie slightly above the scaling curve. The slope 
of the curve, that measures the effective resistivity ex- 
ponent, is correctly reproduced by the theory. We can 
now understand that the deviation of the resistivity ex- 
ponent (cf. Fig. ||a) with respect to its value at the QCP 
(1.5) is due the small but finite value of r. The effective 
exponent only reaches 1.5 for i/A — > co, i.e. for r=0. 
We ascribe the excess amplitude for H = 0, represented 
by the vertical shift, to additional scattering processes 
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dues to disorder. We can compare A, the amplitude of 
the T 2 term in the resistivity in the FL region, with the 
theoretical prediction, A ex l/^/A . The inset in Fig. |l| 
shows the ficld dependence of A as detcrmined from the 
initial slope dp/dT 2 of the resistivity and the theoretical 
prediction. There is good agreement. 
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FIG. 4. Scaling plots of the resistivity along the a-axis. 
The solid line is the theoretical scaling curve. 

The susceptibility can also be calculated and compared 
with the data. The solid line in Fig. (^) is the theoretical 
result for H = 1 kG. The data (and the theoretical curve) 
deviate from a pure T 3 / 4 law as T -> 0. This is due to the 
finite value of r which results in normal FL behavior be- 
low a crossover temperature T FL ~ T A . T FL increases 
with field and can be estimated as sa 0.25 K for H=\ 
T. The crossover to T 2 behavior in x(H — IT) can be 
seen in Fig. (||). The low- field results can be described by 
adding to the dynamical mean-field result an additional 
diverging contribution. The presence of a paramagnetic 
phase giving rise to a T~ 1 divergence of x was suggested 
in Ref. pS} ]. However, to suppress such a contribution 
at 3 K by a field as small as 1 kG one would need that 
the impurities carry huge moments (> 3 hb)- Further- 
more, this hypothesis would not explain the divergence 
of 7 jl6| . The inset in Fig. || shows a fit of the ac data for 
H=l G to the expression X (T) = xmf(T) + clT- - s . This 
can be interpreted in terms of the Griffiths-phase model 
that prediets a power-law divergence 8\ oc T~ 1+A 
with an exponent A that vanishes at the QCP. The value 
A = 0.2 that comes out of our analysis is consistent with 
this picture for a system dose to the QCP. We also found 
that the 7 data Q can be accurately described by the 
analogous expression j(T) — Jmf(T) + a'T~°- 8 where 
7mf(T) oc 1 - bVT is the DMFT prediction @ for C/T. 
The equality of the exponents describing the divergent 
parts of x an( i C/T is one of the predictions of the Gri- 
fith's phase model . 

In summary, we have shown that the NFL properties of 
Ce(Ruo.5Rho.5)2SÌ2 are determined by disorder and prox- 
imity to a QCP. The effeets of the two mechanisms can be 



disentangled by applying a small magnetic field. At low 
fields, disorder effeets dominate. Above lkG, however, 
the temperature and field dependence of the susceptibil- 
ity are well described by the dynamical mean-field theory 
of the spin glass QCP. The T- and i?-dependent resistiv- 
ity is found to obey a universal scaling law. 
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